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We study fermionic zero modes in the domain wall background. The fermions have Dirac and left- 
and right-handed Majorana mass terms. The source of the Dirac mass term is the coupling to the 
scalar field $. The source of the Majorana mass terms could also be the coupling to the scalar 
field $ or the vacuum expectation value of some other field acquired in a phase transition well 
above the phase transition of the field We derive the fermionic equations of motion and find the 
necessary and sufficient conditions for a zero mode to exist. We also find the solutions numerically. 
In the absence of the Majorana mass terms, the equations are solvable analytically. In the case of 
massless fermions a zero energy solution exists and we show that although this mode is not discretely 
normalizable it is Dirac delta function normalizable and should be viewed as part of a continuum 
spectrum rather than as an isolated zero mode. 



I. INTRODUCTION 

The vacuum structure in theories with spontaneous 
symmetry breaking is very rich. There exists topolog- 
ically stabile configurations of gauge and Higgs fields 
known as monopoles, strings and domain walls. Also, 
classical configurations which are not topologically stable 
- non-topological solitons — could exist and be stable 
because of dynamical reasons. Such objects may play an 
important role in the evolution of our universe. 

Extensive research has been done on the interaction 
of fermions with solitons and many new and interesting 
phenomena have been discovered. For example fermionic 
zero modes on strings are responsible for string super- 
conductivity [Q . Quark and leptonic zero modes have an 
important effect on the stability of the non-topological 
electroweak strings ^-Q] . Some work has been done also 
on fermionic zero modes on domain walls, but only Dirac 
fermions have been considered flqjq] . 

Recent experimental results [0 strongly suggest that 
neutrinos, although very light, have a finite mass. Being 
neutral, neutrinos can have Majorana masses in addition 
to the usual Dirac mass terms. This paper study zero 
modes of a particle with both Dirac and Majorana mass 
terms in the background of a domain wall. 

In section 2 we present the Lagrangian of a theory 
where the coupling to the real scalar field $ is the source 
of all mass terms. We derive the equations of motion. 
In section 3 we solve the equations analytically in two 
asymptotic regimes, near the origin and far from the ori- 
gin, and find the necessary and sufficient condition for the 
zero mode to exist. We also solve the equations numeri- 
cally. If the Majorana mass terms are absent it is possible 
to solve the equations analytically and our result agrees 
with the results in literature We discuss the zero 

energy solutions in the case of massless fermions. In sec- 
tion 4 we repeat the procedure of section 3 for the case 
when one or both Majorana masses are spatially homo- 
geneous. The phenomenological significance of obtained 
results is discussed in section 5. 



II. LAGRANGIAN, ANSATZ AND EQUATIONS 
OF MOTION 

We consider a Lagrangian in 3 + 1 dimensions involving 
the real scalar field $ and two chiral spinor fields ip L and 
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where A, rj, k\, ki and k 3 are real and positive (if ipi, and 
i\)r are four component Dirac spinors we can absorb all 
phases of k\, &2 and k 3 into phases of spinor components). 
%jj c = Ctp T , where C is the charge conjugation matrix, is 
introduced in (Q) in order to have the most general mass 
terms for fermions ||. The scalar potential in ([!]) has 
two minima, & m in = ±?7> and in the absence of fermions 
supports the domain wall solution: 



4>o = 77 tanh( 
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for a domain wall located in the xz plane. The resulting 
fermionic equations of motion in a domain wall back- 
ground are: 

i^d^L = h$ ip R + k 3 $ ip c L 
vfdy^R = fci^oV-i + k 2 ^oVJ c R (3) 
The representation of the Dirac matrices we will use 
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where r l , i — 1,2,3 are the Pauli matrices. In 
this representation, a four-component Dirac fermion has 
left and right-handed components of the form: t/r[ = 
(a, /?, -a, -0) and V# = (7, $, 7; 

For a domain wall located in the xz plane, <E>o is a 
function of y only and we look for a solution to (||) of the 
form: 



a = a {y)e iwt - ik * x - ik * z 

7 = c (y)e iu)t ~ ik * x ~ ik * z 
5 = d(y)e iut - ik " x - ik ** 



(5) 



a(y), b(y), c(y) and d(y) can be taken to be real func- 
tions of y. We are interested in zero energy solutions 
to equation (Q) such that all spinor components fall off 
exponentially outside the string core (large y) and are 
well-behaved (nonsingular) at the origin (small y). 

After setting to = k x — k z = 0, we get two identical 
sets of equations: 

a' = /ci$ ff + fc 3 $ a (6) 
d' = /ci$ a + k 2 <fr d 

b' = /ci<J> c+ k 3 <S> b (7) 
c' = k 1 $ b + k 2 $ c 

III. SOLUTIONS 

Let us first analize the system (0) for y > (y < 
is analogous because of the symmetry under y —y)- 
No analytic solution to the system could be found, but 
we can learn something about the structure of solutions 
by looking at their asymptotic behavior. Neglecting a 
domain wall back-reaction to fermions, <E>o — * V at large 
y. Assuming the asymptotic behavior a ~ aoe sy and 
d ~ doe ty (do, do j s an d t are arbitrary real numbers) we 
get following conditions: 

s = t, s 2 - S7j(k 2 +k s ) + k 2 k 3 r, 2 - klrf = (8) 

The solutions to this second order polynomial equation 
for s are 
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(*a + k 3 ) ± (k 2 - k 3 y + 4kj 



(9) 



Both solutions for s are real. s+ is positive, and is giving 
rise to exponentially growing mode. s_ is negative if and 
only if k 2 k 3 < kf, giving rise to exponentially decaying 
mode. Only exponentially decaying modes are physically 
acceptable. So, we conclude that the necessary condition 
for the zero mode to exist on a domain wall is k 2 k 3 < 
k\. We should note that this result is very similar to 



the conjecture in || about the existence of neutrino zero 
modes on electroweak strings. 

The analysis for y < is analogous. For y < we 



have a ~ a$e 



~ s \y\ 



so only positive s leads to a physically 



acceptable result. But we also have 77 * — 77 in (|9|) and 
the necessary condition for the existence of a zero mode 
stays the same. 

If we want to match the exponentially decaying large y 
solutions in a given pair of equations to the solutions at 
the origin, we must know how many solutions for small y 
are well-behaved. Neglecting a domain wall back reaction 
to fermions, $0 ~ foV + /2J/ 3 + ... for small y. Assuming 
that a ~ aoy s , d ~ dgy* (ao and bo are arbitrary real 
numbers while s and t are nonnegative real numbers) for 
small y, we find that the leading orders are s = 0, t = 0. 
System (||) is invariant under parity (y — > —y), so all 
corrections to the solution which are parity violating are 
zero. We can write the general well-behaved solution in 
the form: 



d y° 



a-2y 
d2V 2 



a>4D 
d 4 y 4 
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Substituting (|l(]) into (§) and equating coefficients of 
the same order in y we see that only two coefficients 
are independent (say ao and do); while all others are 
functions of the first two). This means that there are 
two linearly independent well-behaved solutions. For the 
system of two linear first order differential equations we 
should have two solutions in total, so we conclude that 
both solutions are well-behaved. 

Obviously, all the arguments given for y > hold for 
y < too, so we will confine our following discussion 
to positive y. Each of the two well-behaved solutions 
at the origin matches to a unique linear combination of 
both solutions at large y. If both solutions at large y are 
bad (exponentially growing), then there is no solution 
which is well-behaved everywhere. If one of them is good 
(exponentially decaying) there is always one well-behaved 
solution everywhere — good solution at large y which 
matches to a linear combination of two good solutions at 
the origin. In this case, the necessary condition for the 
zero mode to exist, k 2 k 3 < k 2 , is also sufficient. Indeed, 
numerically solving the system (^|), we find that there is 
one well-behaved solution (Figure (|l|)). We took k\r\ = 2, 
k 2 i] = 1 and k 3 r\ = 1 which satisfies the necessary and 
sufficient condition for a zero mode to exist. We also set 
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FIG. 1. Numerical solution of the system (|6|) with k\r] = 2, 
k,2ri — 1, kzr) — 1 and = 1, showing a(y) (solid) and d{y) 
(dashed). The domain wall is located at y = 0. The solution 
is symmetric with respect to y — > — y. 



The system of equations (ffl) have the same structure 
as the system (||) and the previous analysis applies also 
to®. " 

If the fcrmions have only Dirac mass terms, i.e. k\ =/= 
0, k2 = 0, k 3 = 0, it is possible to solve the system 
(^J) analytically. In this case we get the well-behaved 
solution: 



a(y) = A 



cosh(^) 



2fei 



(11) 



cosh( 



where A is a normalization constant. This result is in 
agreement with the results in ||||. Solution sug- 
gests that massless fermions, i.e. fermions with k\ — 0, 
k% = 0, ks = 0, do not have normalizable zero modes. In 
this case, the fermion wave function is a constant and the 
normalization integral diverges. However, as was shown 
in H for the case of massless neutrino zero modes on 
electroweak strings, a correct interpretation of massless 
fermion zero modes is within the framework of the con- 
tinuum. Keeping to in (||) we get equations for massless 
fermions: 



Lua + b' = 
u)b + a' = 
-luc + d! = 
ujd + d = 



(12) 



The solutions are plane waves, a(y), b(y), c(y), d(y) ~ 
e ±KOy , which are Dirac delta- function normalizable. If 



we look at the massless fermion zero mode (us — ► limit 
of a plane wave solution) as an isolated zero mode, it 
is not normalizable. But, this state is actually part of 
a continuum spectrum of the theory and is Dirac delta 
function normalizable. 



IV. CONSTANT MAJORANA MASS TERMS 

Probably the most interesting case is when the source 
of the Majorana masses Ml and Mr is not the coupling 
with the field The Majorana mass terms can be taken 
to be spatially homogeneous and presumed to arise from 
the vacuum expectation value of some field acquired in 
a phase transition above the phase transition of the field 
In this case we set k2$ = Mr and k^ = Ml in 
Lagrangian (|l|). 

We can carry out a similar analysis to the one done in 
section 3. The equations of motion now are: 



a' = fci$ d + M L a 
d' = fci$ a + M R d 

b' = fci$ c + M L b 
c = fci$ fr + M R c 



(13) 



(14) 



An analysis of the large y behavior for y > is analogous 
to the previous case. Substituting fc2?7 = Mr, k^rj = 
Ml in (|^) and demanding negative solutions for s we 
find the necessary condition for a zero mode to exist: 
MlMr < Mp, where we just formally write kirj = Mb- 
For y < 0, we demand positive solutions for s and see 
that if MlMr < Mp there is one decaying mode. 

An analysis of the behavior near the origin is slightly 
different than before since the constant Ml and Mr 
terms break the y — > — y symmetry of the equations (|l3| ) 
and parity breaking corrections to the leading order in a 
solution are also present (odd powers in y). Thus, the 
general well-behaved solution to the system ( |l3|) can be 
written as: 



aoy 

d y° 



aiy 

■diy- 



d 2 y 2 



0-3V 

d 3 y 3 



(15) 



Substituting @ into @ we see that two coefficients are 
independent (say ao and do) which means that there are 
two linearly independent well-behaved solutions. This 
conclusion is valid for both positive and negative y. The 
matching of well-behaved solutions at large and small y is 
similar to system (^). However, the matching of a given 
pair of solutions for a(y) and d(y) for y > and y < 
at y — is nontrivial. There is no symmetry y — ► — y 
¥= d(^o-) ( Fi S ure (§))• Numerical analy- 
only if M L = Mr. We 
conclude that if MlMr < Mf,, where Ml = Mr, there is 
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one solution which is well-behaved everywhere. Numer- 
ically solving the system dl3| ) setting kit] = 2, Mr = 1, 

Ml = 1 and = 1, we find one well-behaved solution 
(Figure (§)). 
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FIG. 2. Numerical solution of the system ( |l3| ) with 
km = 2, M R = 1, Mi = 1.5 and ^f 2 = 1, showing a(y) 
(solid) and d(y) (dashed). We see that Jjg|=g±} / 
and the solution is not continuous. 
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FIG. 3. Numerical solution of the system ( |l3| ) with 

fcuy = 2, Mi? = 1, Ml = 1 and = 1, showing a(y) (solid) 
and d(y) (dashed). The domain wall is located at y — 0. The 
solution is not symmetric with respect to y — * —y. 

Let us mention that a similar analysis can be done 
in the case when just one Majorana mass term is ho- 
mogeneous, while the second one still comes from a 
coupling to the field $. Setting fc 3 $ = M L in (Q) 
and repeating the procedure from section 3 we find the 



same form of the necessary condition for a zero mode 
to exists M l Mr < Ml), where we just formally write 
kit] = Mjj and &2?7 = Mr. However, unlike the previ- 
ous case, we can not set M L = k 2 $o everywhere. Thus, 
% y= °t\ + a Ji y=0 n ~\ for a11 values of m l and there is no 

d(y=0+) ' d(y=0-) u 

matching between the y > and y < solutions for a(y) 
and d(y). In this case there are no well behaved solutions. 

V. PHENOMENOLOGY 

By now, we have not addressed the question about the 
phenomenological validity of the general fermionic mass 
terms. The Lagrangian (0) could be, in principle, part of 
some larger theory like the standard model (and various 
extensions of it) or GUT theories. The specific charges 
assigned to the fermionic field ifi would then depend on 
a specific symmetry group of the model. Mass terms, 
present in a Lagrangian have to be gauge singlets, thus 
restricting their possible forms. An illustrative exam- 
ple, phenomenologically valid standard model neutrinos, 
beside the Dirac mass terms, can only have the right- 
handed Majorana mass terms. In our notation, this cor- 
responds to the case of k 3 = 0. If the coupling to the 
scalar field $ gives rise to the right-handed Majorana 
mass we analize the results of section 3. Setting &3 = 
in (||) we see that the necessary condition k 2 k^ < kf, 
i.e. < fcf, is always satisfied. Following the rest of the 
analysis of section 3, we conclude that in this case there 
is always one zero mode. However, in this case it would 
be difficult, without a fine tunning, to explain a large 
difference between the Dirac and right-handed Majorana 
masses because the same phase transition gives the scale 
for both masses. 

If the right-handed Majorana mass is spatially homoge- 
neous, then we apply the analysis of the section 4 (which 
is different from the analysis of the section 3 because 
there is no symmetry y — > —y). In the absence of the 
left-handed Majorana mass, although the necessary con- 
dition is satisfied, due to the fact that = Ml ^ Mr, 
the solutions for positive and negative y do not match 
at y = 0. In this case, there are no well behaved zero 
modes. 

VI. CONCLUSION 

We studied zero energy solutions of the Dirac equa- 
tion in the background of a domain wall located in the 
xz plane. We first considered the case when the vacuum 
expectation value of the scalar field $ was the source 
of the Dirac and Majorana mass terms. In the general 
case, when both the Dirac and Majorana mass terms are 
present we solved the equations of motion analytically in 
two asymptotic regimes, large and small y and found the 
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necessary and sufficient condition for a zero mode to ex- 
ist. This condition is very similar to the one conjectured 
in J! in the case of neutrino zero modes on electroweak 
strings. We also solved the equations numerically. 

If the Majorana mass terms are absent, it is possible to 
solve the equations analytically. The result agrees with 

§§■ 

If the fermions are massless there exists a zero-energy 
solution which is not discretely normalizablc. Solving 
the non-zero energy equations we showed that this state 
is actually part of the continuum spectrum and is Dirac 
delta function normalizable. 

We also considered the cases when one or both Majo- 
rana terms are spatially homogeneous. If both of them 
are spatially homogeneous and if Ml Mr < M|>, with 
Ml = Mr, there is one well behaved solution. If just 
one is spatially homogeneous there are no well behaved 
solutions. The phenomenological significance of obtained 
solutions was discussed. 
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